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Abstract 

We give blueprints of the Schwinger-DeWitt technique for the covariant curvature ex- 
pansion of quantum effective action for the DGP type models in curved spacetime. 

1 Introduction 

Modified theories of gravity in the form of braneworld models can in principle account for 
the phenomenon of dark energy as well as for nontrivial compactifications of multi-dimensional 
string models. It becomes increasingly more obvious that one should include in such models 
the analysis of quantum effects beyond the tree-level approximation pQ. This is the only way 
to reach an ultimate conclusion on the resolution of such problems as the presence of ghosts [2 J 
and low strong-coupling scale [3]. Quantum effects in brane models are also important for the 
stabilization of extra dimensions [3J, fixing the cross-over scale in the Brans-Dicke modification 
of the DGP model [5] and in the recently suggested mechanism of the cosmological acceleration 
generated by the four-dimensional conformal anomaly [BJ. 

A general framework for treating quantum effective actions in brane models (or, more gen- 
erally, models with timelike and spacelike boundaries) was recently suggested in [7} 151 |9"1 [TU]. 
The main peculiarity of these models is that due to quantum field fluctuations on the branes 
the field propagator is subject to generalized Neumann boundary conditions involving normal 
and tangential derivatives on the brane/boundary surfaces. This presents both technical and 
conceptual difficulties, because such boundary conditions are much harder to handle than the 
simple Dirichlet ones. The method of [9] provides a systematic reduction of the generalized 
Neumann boundary conditions to Dirichlet conditions. As a byproduct it disentangles from the 
quantum effective action the contribution of the surface modes mediating the brane-to-brane 
propagation, which play a very important role in the zero-mode localization mechanism of the 
Randall-Sundrum type [TT]. The purpose of this work is to make the next step — to extend a 
well-known Schwinger-DeWitt technique |12|, [13l [Pi} [T5] to the calculation of this contribution in 
the DGP model in a weakly curved spacetime in the form of the covariant curvature expansion. 

Briefly the method of [9] looks as follows. The action of a (free field) brane model generally 
contains the bulk and the brane parts, 

S[4>\ = \J d d+1 X G 1 / 2 0(A)F(Vx) 4>{X) + \J d d x g x l\{x) «(V X ) <p(x) , (1) 

where the (d+ l)-dimensional bulk and the (i-dimensional brane coordinates are labeled re- 
spectively by X = X A and x = x^, and the boundary values of bulk fields (j>(X) on the 
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brane/boundary b = SB are denoted by ip(x), 



cP(X) =<p(x), (2) 

b 

G and g are the determinants of the bulk Gab and metrics respectively. 

The kernel of the bulk Lagrangian is given by the second order differential operator F(Vx), 
whose covariant derivatives Vx are integrated by parts in such a way that they form bilinear 
combinations of first order derivatives acting on two different fields. Integration by parts in the 
bulk gives nontrivial surface terms on the brane/boundary. In particular, this operation results 
in the Wronskian relation for generic test functions <f>i 2{X), 

d d+1 X G 1 / 2 ( 0! F(V x)4>2 - 4>iF(V x ) <h) = ~ J d W /2 (faWfa-fa W<h) • (3) 

Arrows everywhere here indicate the direction of action of derivatives either on 0i or 02 • 

The brane part of the action contains as a kernel some local operator k(V), V = V x . Its 
order in derivatives depends on the model in question. In the Randall- Sundrum model [llj . 
for example, it is for certain gauges just an ultralocal multiplication operator generated by the 
tension term on the brane. In the Dvali-Gabadadze-Porrati (DGP) model [16] this is a second 
order operator induced by the brane Einstein term on the brane, re(V) ~ VV/m, where m is 
the DGP scale which is of the order of magnitude of the horizon scale, being responsible for 
the cosmological acceleration [T7]. In the context of the Born-Infeld action in D-brane string 
theory with vector gauge fields, k(V) is a first-order operator |18j . 

In all these cases the variational procedure for the action (pQ) with dynamical (not fixed) 
fields on the boundary <p(x) naturally leads to generalized Neumann boundary conditions of 
the form 

'w(v x ) + K(yj)<i> =o, (4) 

/ b 

which uniquely specify the propagator of quantum fields and, therefore, a complete Feynman 
diagrammatic technique for the system in question. The method of [9] allows one to sys- 
tematically reduce this diagrammatic technique to the one subject to the Dirichlet boundary 
conditions <j>\ b = 0. The main additional ingredient of this reduction procedure is the brane 
operator F hTa,nc (x,x') which is constructed from the Dirichlet Green's function Gd(X,X') of 
the operator F(V) in the bulk, 

F b ^(x,x') = - W(V X )G D (X,X')W(V X/ ) +k(V)5(x,x') . (5) 

X=e(x), X =e[x ) 

This expression expresses the fact that the kernel of the Dirichlet Green's function is being 
acted upon both arguments by the Wronskian operators with a subsequent restriction to the 
brane, with X = e{x) denoting the brane embedding function. 

As shown in [9], this operator determines the brane-to-brane propagation of the physical 
modes in the system with the classical action ([T|) (its inverse is the brane-to-brane propagator) 
and additively contributes to its full one-loop effective action according to 

A-loop = \ Tr!J +1) InF = i Trg +1) hxF + \ Tr^ lnF brane , (6) 

where Tr^ 4 ^ denotes functional traces of the bulk theory subject to Dirichlet and Neu- 
mann boundary conditions, respectively, while Tr^ is a functional trace in the boundary d- 
dimensional theory. The full quantum effective action of this model is obviously given by the 
functional determinant of the operator F(Vx) subject to the generalized Neumann boundary 
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conditions ©, and the above equation reduces its calculation to that of the Dirichlet boundary 
conditions plus the contribution of the brane-to-brane propagation. 

Here we apply ([6]) to a simple model of a scalar field which mimics in particular the properties 
of the brane-induced gravity models and the DGP model [16]. This is the (d+l)-dimensional 
massive scalar field <j)(X) = 4>(x, y) with mass M living in the curved half-space y > with the 
additional d-dimensional kinetic term for (p{x) = (j)(x, 0) localized at the brane (boundary) at 

y = o, 

S[<j>} = \ f ^^((V^Ilf + MVtXjj+i j d<W /2 (V^(*)) 2 . (7) 

Here and in what follows we work in a Euclidean (positive-signature) spacetime. Therefore, 
this action corresponds to the following choice of F(Vx) m terms of (d+ l)-dimensional and 
d-dimensional covariant D'Alembertians (Laplacians) 

F(V X ) = M 2 -D( d+1 ) = M 2 -G AB V A V B ,- (8) 

In the normal Gaussian coordinates its Wronskian operator is given by W = —d y — the normal 
derivative with respect to outward-pointing normal to the brane, and the boundary operator 
k(V) equals 

«(v) = ~n, n = nW=^v M v„ (9) 

where the dimensional parameter m mimics the role of the DGP scale [16] . Thus, the generalized 
Neumann boundary conditions in this model involve second-order derivatives tangential to the 
brane, 



0, (10) 



cf. dH) with W = -d y and k = -D/2m. 

As was shown [lOj . the flat space brane-to-brane operator for such a model has the form of 
the pseudodifferential operator with the flat-space □, 

F branc (V) = — (-□ + 2my / M 2 - □). (11) 
2m 

In the massless case of the DGP model [16] . M = 0, this operator is known to mediate the 
gravitational interaction on the brane, interpolating between the four-dimensional Newtonian 
law at intermediate distances and the five-dimensional law at the horizon scale ~ 1/m [3\. 

Here we generalize this construction to a curved spacetime and expand the brane-to-brane 
operator and its effective action in covariant curvature series. This is the expansion in powers 
of the bulk curvature symbolically denoted below as R, extrinsic curvature of the brane denoted 
by K and their covariant derivatives — all taken at the location of the brane. The expansion 
starts with the approximation (jlip based on the full covariant d'Alembertian on the brane. We 
present a systematic technique of calculating curvature corrections to (jlip and rewrite their 
nonlocal operator coefficients — functions of the covariant □ — in the form of the generalized 
(weighted) proper time representation. 

The success of the conventional Schwinger-DeWitt method is based on the fact that the 
one- loop effective action of the operator, say M 2 — □, has a proper time representation 

^Trln (M 2 -n) = -\ f° ^ e~ sM " Tr e sU . (12) 



3 



In view of the well-known small time expansion for the heat kernel [12\ [T3] , 

^ oo 

e sU 5{x,x>) = ——-jj^D 1 / 2 (x, x') e -<r(. x > x ')/ 2a ^ s n a n (x,x% (13) 

(<t(x,x') is the geodetic world function, D(x,x') is the associated Van Vleck determinant and 
a n (x,x') are the Schwinger-DeWitt or Gilkey-Seely coefficients) the curvature expansion even- 
tually reduces to the calculation of the coincidence limits of a n (x,x') and a trivial proper time 
integration resulting in the inverse mass expansion 

^ ' 71=0 

As we will show below, the calculation of the brane effective action differs from the conven- 
tional Schwinger-DeWitt case in that the proper time integral (|12p contains in the integrand a 
certain extra weight function w(s) and instead of just Tre s ^ one has to calculate the trace of 
the heat kernel acted upon by a certain local differential operator Tr (VF(V)e sC ) . This again 
reduces to the calculation of the coincidence limits — this time of the multiple covariant deriva- 
tives of a n (x, x'), V m ...V fln a n (x, x')\ x / =x — the task easily doable within a conventional DeWitt 
recurrence procedure for a n (x,x'). 

2 Perturbation theory for the bulk Green's function and brane 
effective action 

In normal Gaussian coordinates the covariant bulk d'Alembertian decomposes as □ x = 
dy + D(y) + where ellipses denote depending on spin terms at most linear in derivatives^ and 

□ (y) is a covariant d'Alembertian on the slice of constant coordinate y. Therefore the full bulk 
operator takes the form 

F(V) = M 2 - D { x +1) + P{X) = M 2 — H — dy — V(X \ d y ,V) = F° — V, (15) 

in which all nontrivial y-dependence is isolated as a perturbation term V(X \ d y ,V) = V(y, d y ) 
— a first-order differential operator in y, proportional to the extrinsic and bulk curvatures, and 
of second order in brane derivatives V which we do not explicitly indicate here by assuming that 
they are encoded in the operator structure of V(y, d y ). In particular, it includes the difference 

□ (0) — □(?/) = □ — 0(y) expandable in Taylor series in y. 

The kernel of the bulk Green's function can formally be written as a y-dependent nonlocal 
operator acting on the (i-dimensional brane — some non-polynomial function of the brane 
covariant derivative 

G D (X,X') = G D (y,y'\V)5(x,x'). (16) 
The perturbation expansion for Go(y, y'\ V) is usual 

oo 

G D = G% + G%VG% + ... = G° D £ (G% V) n , (17) 

n=0 

where G® D is the propagator for operator F° obeying Dirichlet boundary conditions and the 
composition law includes the integration over the bulk coordinates, like for example in the first 
subleading term 

/•oo 

G D VG° D (y,y')= / dy" G° D (y, y")V(y", d y ii)G° D (y" , y). (18) 
Jo 



1 This term for a general spin structurally has the form KVx + K 2 + (VA") + R where K is the extrinsic 
curvature of y — const slices and R is the bulk curvature. 
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The lowest order Green's function in the half-space of the DGP model setting — the Green's 
function of F° = M 2 — □ — d 2 subject to Dirichlet conditions on the brane y = and at infinity 
— reads as follows 

-| y~y'\y r W^U _ p -(y+y') y r W^D 

G° D (y,y') = - ; . (19) 

m ' 2VM 2 - □ V ' 

We want to stress that here we assume the exact (curved) d-dimensional d'Alembertian □ 
depending on the induced metric of the brane g tiv {x). This means that in the lowest order 
approximation the underlying spacetime is not flat, but rather has a nontrivial but constant in 
y metric of constant y slices. Correspondingly in the zeroth order we have 



[WG D (y,y')W]° Q = d y G° D (y, y') d y = -VM 2 -D. (20) 

y y y=y'=Q 

The perturbation of the bulk operator can be expanded in Taylor series in y, so that it reads 



OC 



V(y,d y ) = ^2y k V k (d y ), (21) 

fc=0 

where Vt (d y ) = Vt{d y \SJ) is a set of y-independent local d- dimensional covariant operators of 
second order in V x and first order in d y . 

On substitution of (|20"1) and (|2T1) into (fTTl) exactly calculable integrals over y result in a 
nonlocal series in inverse powers of \JM 2 — □, and the perturbation expansion takes the form 



[WG D {y,y')w] y=yf=Q = - VM 2 ~ □ + £ U k ( V) — _ g)fc/2 , (22) 

where C/fc(V) is a set of local covariant differential operators acting on the brane^- The di- 
mensionality of each f7^(V) is the inverse length to the power k+1, which is composed of the 
dimensionalities of bulk and extrinsic curvatures and covariant derivatives all taken on the brane 
at y = 0. 

With k(V) = — D/2m the brane-to-brane operator reads 

oo _. 

2mF b r anc (v) = _ n + 2my / M2 _ D _ 2m ^ ^ (V ) — _ Q)fc/2 , (23) 

Then we consider the perturbation series for the functional trace of its logarithm in powers of 
the full [/^-series. After reexpansion in powers of two sets of nonlocal propag ators 1/VM 2 - □ 
and l/(— □ + 2m \/ M 2 — □) the brane effective action finally takes the form 



TrlnF brane = -Trln ( - □ + 2mV ' M 2 - □ 

+ ^ W »W (M > - U)^ (-□ + Wtf - □)■ P4) 

with a new set of local covariant differential operators Wki(V) acting on the brane. The dimen- 
sionality of Wjw(V) is k + 21 in units of inverse length. One should also remember that each 
power of l/(— □ + 2m\/M 2 — □) is accompanied by one power of m in the numerator, so that 
structurally 

W W (V) ~ m l X7 a R b K c , (25) 

where the integer overall powers of the covariant derivatives V, bulk curvatures R and extrinsic 
curvatures K are constrained by the relation a + 2b + c = k + I. 



2 Strictly speaking each fc-th order of this series arises in the form of the following nonlocal chain of square root 
"propagators", (M2 _ 1 D) i 1 /2 Ui (M2 _ 1 D) i 2 /2 U 2 ...U p -i ^ lp/2 , h + h + ■■■ + l p = k, with differential operators 
Ui as its vertices, but all these propagators can be systematically commuted either to the uppermost right or left 
by the price of extra commutator terms of the same structure. 
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3 Generalized proper time method 



Our goal now is to find the proper time representation of nonlocal operators in Eq.(|24 
in the form of the exponentiated □. A systematic way to do this consists in the following 
factorization of the brane-to-brane operator as 



2mF brane (V) = -□ + 2m\/M 2 - □ = (\/ M 2 - □ - m + )(V ' M 2 - □ 



Here the masses m± are the roots of the relevant quadratic equation, x 2 + 2mx — M 2 
x = VM 2 - □, 



m± 



-m 



± \JM 2 + m 2 , m_ < —M < < m + < M, 



(26) 
= 0, 

(27) 



which determine the poles of the propagator of i* 1 Q ranc (V) in spacetime with the Lorentzian 
signatur^l 



□_ 



M 
\M 2 



m\ > 0, 



m 2 _\e 3in < 0. 



(29) 
(30) 



The factorization (|26|) allows one to rewrite the l-th. power of the brane-to-brane propagator 
in (El as 



1 



1 



1 



(-□ + 2m\/M 2 - □ ) l {\/W 2 



□ 



mi 



(VM 2 - □ 



rri- 



(31) 



and then decompose the resulting fraction into the sum of simple fractions for which one has 
explicit proper time representations. These representations begin with the following relations 

m 



i 



i 



(M 2 - ufi 2 r(k/2) 



ds s k/2-l e s(n-M*) 



1 



\/M 2 - U-m 



oo 

J ds e s ( Q - M2 ) (^j^ + mw{-myfs) S j , m < M, 



VM 2 - □ (VM 2 - □ 



m 



wl—m^/s), m < M, 



(32) 
(33) 
(34) 



which generate (by differentiating with respect to m, □, M and linear recombining the results) 
the list of fractions with all possible powers of the factors VM 2 - □ and yjM 2 - □ — m± m 
their denominators. Here the weight function w{s) is given in terms of the error function [19] 
erf(x) = Jq dy exp(— y 2 ) and has the following ultraviolet and infrared asymptotics 



w(x) 



e x2 [l 



erf (x) 



1. 

l 



2e a 



x 
x 

X 



0, 

+oo, 
— oo 



(35) 



3 The pole at □_ is formally tachyonic, but it is located on the unphysical sheet of the Riemann surface for 
the propagator in the complex plane of □ [20] (which is indicated in (|3(jp by the nontrivial phase) . Moreover, its 
residue is identically vanishing in view of m_ < 0. Therefore this pole does not correspond to a real particle. For 
M^0 only □+ gives rise to a particle with the decreasing mass as M — + 0, Af 2 — rn\ — * 0, which also disappears 
in the DGP limit because the pole residue also vanishes at M = 0. In this limit only the continuum spectrum of 
massive intermediate states survives forming the spectral representation for the DGP propagator [21] 



-□ + 2m V c n 



4m f c 
7T Jo 



1 



+ 4m 2 p 2 — □ ' 



(28) 
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The last two proper time integrals above are defined only for m < M (for any negative m 
and for m < M if m is positive), because in view of these asymptotics they are convergent at 
infinity only in this range. Interestingly, the forbidden domain corresponds to the real tachyon, 
because for m > M > the pole □ = M 2 — m 2 belongs to the physical sheet of the propagator, 
and its residue is nonvanishing. 

From (|33j) - ([34j ) it immediately follows that the zeroth order brane-to-brane propagator and 
its one-loop functional determinant read 



J = f dse s (D-M 2 ) m+w(-m +y /s) - m-w(-m-y/s) 

J m+ — m- 



-□ + 2mVM 2 - □ 

o 



OC 



Tr In (- □ + 2m\/M 2 - □ 

= -\tt J ^-e s ^ n - M ^ (w(-m + ^)+w(-m.^, (37) 
o 

For the DGP model case with M 2 = and m+ = 0, m_ = —2 m these representations 
simplify to the equations derived in [10] 



-□ + 2m v /z O 



J dse sD w(2m^/s), (38) 
o 

OO _ 

Tr In - □ + 2mv /Z S = -Tr / — e s u i —. (39) 



s 2 



o 



The interpretation of the weight contribution here is very transparent. It interpolates between 
the ultraviolet and infrared domains where the brane operator and its logarithm have quali- 
tatively different behaviors. In the domain of a small proper time m^/s <C 1 it approximates 
the brane operator by a large — □ S> m 2 (hence the overall weight (1 + w)/2 ~ 1), whereas in 
the infrared domain m^fs S> 1 it approximates the operator by 2my/— □ (hence the weight is 
(l + io)/2~ 1/2 corresponding to ln v /= n= (1/2) ln(-D) ). 

By decomposing the nonlocal fractions of (|24j) into the sum of simple fractions and using 
the weighted proper time representations (I32p - (|34p and their derivatives with respect to mass 
parameters m± we obtain the following expression 



OO 

i i r ds 



(m 2 - n) fc / 2 (-□ + 2mVM 2 -ay J s 



e 



■ 8(MJ - D) %( S ,m,M), (40) 



o 



with some weight function Wki{s,m, M)@ 

Using (|37p and (|40p we finally obtain for the perturbative expansion (|24p 



- TrlnF brane = -- [ — e~ sM2 w (~ m +v /i ) + w(-m-y/s) sD 
2 2 J s 2 

o 

OO 



+ £ f^-e~ sM2 w kl ( S ,m,M)Tt[w M (V)e sD ]. (41) 

k>0,l>l o 



4 Alternatively this weight function can of course be obtained as a Mellin transform of the function of □ in 
the left hand side, but this simple fraction decomposition method gives a more regular and systematic way to 
achieve the needed goal. 
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This formally solves the problem of constructing the Schwinger-DeWitt expansion for the brane 
effective action, because as it was expected all the remaining calculations reduce to the con- 
ventional calculation of the coincidence limits of the Schwinger-DeWitt coefficients and their 
covariant derivatives in 

ft N V)e 1 = (i/ ^,V 2 f; s n tr ^(v^-^^)a„(x,,')|^. (42) 

^ ' n=0 

Remember that every Wjy(V) is a finite order covariant differential operator with the coefficients 
built of the powers of the bulk curvature, extrinsic curvature of the brane and their covariant 
derivatives. Here lengthening of the derivatives in Wfcj(V) originates from commuting them 
with the the exponential factor exp(— a(x, x')/2s) contained in the kernel of exp(sD), a^(x, x') = 
V*a(x, x'). This of course brings to life world function coincidence limits V w ...V^ p a(x, x')\ x > =x 
also easily calculable by the DeWitt recurrence procedure [El [13] . 

It is important that the expansion (|4ip is efficient for the purpose of obtaining the asymptotic 
1/M-expansion. Indeed, in view of the weight function asymptotics (|35p the w(— m_\/i)-parts 
of the overall wm(s, m, M) (cf. Eq. ([36]) ) with m_ < are suppressed at s — > oo by e~ sM and, 
therefore, generate after the integration over s the needed l/M 2 -series. In the w(— m + yfs)- 
parts the integrand behaves as e~ s ^ M ~ m +\ and generates the 1/(M 2 — m^_) ~ l/2?nM-series 
also appropriate for the M — ► oo limit, though converging slower than the 1/M 2 one. This 
is the expansion in inverse squared masses of the real particle associated with the pole (|29p . 
Unfortunately, the powers of 1/M are accompanied by those of 1/m, which comprises in the 
DGP model the problem of low strong-coupling scale [3] for small DGP crossover scale m. 



4 Conclusions 

This is obvious that the Schwinger-DeWitt technique in brane models is much more compli- 
cated than in models without spacetime boundaries. It does not reduce to a simple bookkeeping 
of local surface terms like the one for simple boundary conditions reviewed in [15j . Neverthe- 
less it looks complete and self-contained, because it provides in a systematic way a manifestly 
covariant calculational procedure for a wide class of boundary conditions including tangential 
derivatives (in fact of any order). On the other hand, the calculational strategy of the above 
type is thus far nothing but a set of blueprints for the Schwinger-DeWitt technique in brane 
models, because there is still a large set of issues and possible generalizations to be resolved in 
concrete problems. 

One important generalization is a physically most interesting limit of the vanishing bulk 
mass M 2 . Local curvature expansion is perfect and nonsingular for nonvanishing M 2 and 
applicable in the range of curvatures and magnitudes of spacetime derivatives (R, K 2 ,VK) <C 
M 2 , VV-R <C M 4 , etc. However, for M 2 — ► it obviously breaks down, because the proper 
time integrals start diverging at the upper limit. These infrared divergences can be avoided by 
a nonlocal curvature expansion of the heat kernel of [22J. Up to the cubic order in curvatures 
this expansion explicitly exists for Tr e sD [23J, but for the structure involving a local differential 
operator TrW(V)e s ^ it still has to be developed. 

Another important generalization is the extension of these calculations to the cases when 
already the lowest order approximation involves a curved spacetime background (i.e. dS or 
AdS bulk geometry, deSitter rather than flat brane, etc.). The success of the above technique 
is obviously based on the exact knowledge of the y-dependence in the lowest order Green's 
function in the bulk and the possibility to perform exactly (or asymptotically for large M 2 ) the 
integration over y. All these generalizations and open issues are currently under study. 

To summarize, we developed a new scheme of calculating quantum effective action for the 
braneworld DGP-type system in curved spacetime. This scheme gives a systematic curvature 
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expansion by means of a manifestly covariant technique. Combined with the method of fixing 
the background covariant gauge for diffeomorphism invariance developed in [HI [23] this gives 
the universal background field method of the Schwinger-DeWitt type in gravitational brane 
systems. 
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